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. Abstract 

CN 

Based on the assumption that an S4 flavor symmetry is embedded into SU(3), a lepton 
^ ' mass matrix model is investigated. A Froggatt-Nielsen type model is assumed, and the flavor 

structures of the masses and mixing are caused by VEVs of SU(2)i,-singlet scalars (j)u and (pd 
which are nonets (8+1) of the SU(3) flavor symmetry, and which are broken into 2+3+3' and 
1 of S4. If we require the invariance under the transformation {(j)'-^'> , cj)'-^^) (— +0^-'^^) 
' (— I 'I for the superpotential of the nonet field 0^^+^), the model leads to a beautiful relation for 

CIh' the charged lepton masses. The observed tribimaximal neutrino mixing is understood by 

pi • assuming two SU(3) singlet right-handed neutrinos i^j^' and an SU(3) triplet scalar x- 
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1 Introduction 

The observed mass spectra and mixings of the fundamental particles will provide promising 



> , 

10 ! clues to unified understanding of the quarks and leptons. Especially, in the lepton sector, the 

, following characteristic features have been observed [I]: 

^ ', (i) The observed charged lepton masses {me,m^,mr) satisfy the relation [21 [3] 
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with remarkable precision; 
(ii) The o 
mixing [3] 
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5-H I (ii) The observed neutrino mixing Ui, is approximately given by the so-called tribimaximal 
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Such characteristic features have not been seen in the quark sector. For example, the mixing 
form (1.2) suggests that the mixing can be described by Clebsh-Gordan-like coefficients, while, 
for the Cabibbo-Kobayashi-Maskawa mixing in the quark sector, such a characteristic feature 
has not been seen, although we have known some relations among the mixing angles and quark 
mass ratios. Therefore, for a start, in the present paper, we investigate the lepton masses and 
mixings. 

In order to understand the relation (1.1), for example, we assume that there are three 
scalars (pi {i = 1,2,3), and the values of the charged lepton masses niei are proportional to the 
square of the vacuum expectation values (VEVs) Vi = {(pi) of the scalars (pi, rriei = kvf (in the 



Ref . [3l El [6] , for instance, a seesaw type model {M^jij = 5ijVi{ME) ^vj has been assumed). We 
define singlet (pa and doublet {(p-jr^cpr]) of a permutation symmetry S3 [7j by 
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from the three objects {(I)i,(p2,4'3): and we consider the following S3 invariant scalar potential 

m mum- 



The minimizing condition of the potential (1.4) leads to the relation 
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The relation (1.5) means 
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Therefore, we can obtain the mass relation (1.1). Here, note that although the scalar potential 
(1.4) is invariant under the S3 symmetry, but it is not a general one of the S3 invariant form. As 
pointed out in Ref. [9j, the scalar potential with a general form cannot lead to the relation (1.5). 
For the derivation of the VEV relation (1.5), it is essential to choose the specific form (1.4) of 
the S3 invariant terms. Similar formulation is also possible for other discrete symmetries A4 |10j 
and S4 (see below). However, in such a symmetry, we still need an additional specific selection 
rule. What is the meaning of such a specific selection? In the present paper, we investigate this 
problem by assuming that the S4 flavor symmetry is embedded into SU(3). 

Recently, a superpotential which leads to the relation (1.5) has proposed by Ma [Hj on the 
basis of a symmetry S(81). Stimulated by the Ma's idea, the author [10] has also investigated 
a similar superpotential on the basis of a symmetry A4. Here, based on an S4 flavor symmetry 
instead of the A4 symmetry, let us review the superpotential W which gives the relation (1.5). 
We denote singlet and doublet of S4 as (pa and (po = {cp-j^ , cprj)'^ , respectively, as well as those 
in S3. In order to write the superpotential for the scalar fields (pa and doublet (pu of S4, we 
put the following phenomenological rule [TU]: the field (pa {a = a,D) to the power nth, {(pa)"' 
(n = 1,2,3), appears always accompanied with the factor 1/n! in the superpotential W. Under 
this phenomenological rule, we can uniquely write the superpotential of (pa and (pn sts 





The potential (1.8) can also lead the relation (1.5). What is the meaning of this phenomenological 
rule? 

On the other hand, we have to consider a mechanism which yields the charged lepton masses 
vrLei oc vf, i.e. the effective Hamiltonian for the charged lepton sector 

Hl^f = [eLi{(t)ifeRi + eL2{(t>2feR2 + eizickzfem] . (1.9) 

We will propose a Froggatt-Nielsen type model [12j, Hl^^ = {IiHfcfxfteji) in Sec. 4. 

Now, let us return the topic of the tribimaximal mixing. From the definition (1.2), we can 
denote the fields ("01 j "02, "03 ) as 
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The observed neutrino mixing (1.2) means that when the mass eigenstates of the charged lep- 
tons are given by the ('i/'i) V'2, V's) basis, the mass eigenstates of the neutrinos are given by the 
{ipn,'4'a,ipn) basis. Therefore, the problem is to find a model where the charged lepton mass 
eigenstates are (ei, 62, 63), while the neutrino mass eigenstates are given by (i/,,, i/g-, i^n) with the 
masse hierarchy < <^ (or <C rn^ < m^). In the present paper, we will investigate 
such a model based on an S4 model. Here, note that the fermions (^/^i, -02) V^s) is a triplet of 
S4, but the basis (■0,^, V'o-, ■07r) is not in any irreducible representations of S4, while the scalar 
(01,02)03) is not irreducible representation of S4, but (07r,0,7) and (p^j are doublet and singlet 
of S4. 

Thus, the characteristic features (1.1) and (1.2) in the lepton sector may be understood from 
the language of S4 (also S3 or A4). However, as seen from the above review, the characteristic 
features (1.1) and (1.2) cannot be understood from the S4 symmetry only. We need some 
additional assumptions. In this paper, we will investigate these problems under an assumption 
that the present S4 symmetry is embedded into an SU(3) symmetry |13j . In the next section, 
the singlet 00- and doublet (0,^, 4>ri) will be understood as members of a nonet scalar [1+8 of 
SU(3)], and the VEV relation (1.5) will be derived by requiring that W^(0) is invariant under a 
Z2 symmetry. 

We know that the three masses in any sectors of quarks and leptons are completly different 
among them. Therefore, if we assume a fiavor symmetry, the symmetry must finally be broken 
completely. Usually, a relation which we derive in the exact symmetry limit is only approximately 
satisfied under the symmetry breaking. Although we derive the VEV relation (1.5) under the 
S4 symmetry, the problem is whether the VEV relation (1.5) which is obtained under the S4 
symmetry is spoiled or not when we introduce such a symmetry breaking. In Sec. 3, we will 
demonstrate that such a symmetry breaking term without spoiling the relation (1.5) is indeed 
possible. 



In Sec. 4, in order to give the charged lepton masses and tribimaximal neutrino mixing, we 
win discuss the effective Hamiltonian by assuming an Froggatt-Nelsen |12] type model. FinaUy, 
Sec. 5 wiU be devoted to the summary and concluding remarks. 

2 VEVs of SU(3) nonet scalars 

The goal in the present section is to obtain the VEV relation (1.6) [i.e. (1.5)]. As seen 
in the previous section, in order to obtain the desirable results (1.5), we need assume an equal 
weight between the doublet and singlet terms of S4. In the present paper, we assume that the 
S4 symmetry is embedded into an SU(3) symmetry. The doublet ((/>7r, (/>7j) and singlet (f)^ of S4 
are embedded in the 6 and (8 + 1) of SU(3) [13]. In the present paper, we assume that the 
doublet {(pTr,4'v) singlet (pa- originate in SU(3) octet and singlet, respectively. The essential 
assumption in the present paper is that the fields (pu and (p^i always appear in the theory with 
the form of the nonet of U(3): 
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and the index f {f = u, d) has been dropped. 

The outline to obtain the superpotential form (1.8) in the present scenairo is as follows: 
The SU(3) invariant superpotential for the nonet fields 4>f {f = u,d) are given by 



WicPf) = ^mfTT{<Pf<t>f) + ^ 



A/Tr(0/0/</./). 



(2.3) 



Since, in the next section, we want to assign chages +1 and — 1 of a Z3 symmetry to the fields 
(pu and 4>d, respectively, we also assign the Z3 charges +1 and —1 to the mass parameters m„ 
and rud in Eq.(2.3), respectively. However, since we do not consider a mass term Tr^cpu'pd), we 
do not consider a mass parameter with the Z3 charge zero. Hereafter, in the present section, for 
convenience, we will drop the index /, since the cross terms between 0„ and (pd do not appear. 



In the superpotential (2.3), although the term Tr 
of S4, the cubic term Tr{(p(p(p) gives 



gives the desirable term (pi^ + + (pi^ + ■ 
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where the terms "• • • " denote terms which include 3 and 3' of the subgroup S4. Therefore, the 
potential (2.3) cannot give the relation (1.5). We must drop the first term in the cubic terms 



(2.4). For this purpose, we introduce a Z2 symmetry, and we assign the Z2 parities —1 and 
+1 (the Z2 charge +1 and 0) for the octet part (f)^^^ and singlet part (f)^^^ of the nonet field (f), 
respectively. The symmetry Z2 breaks U(3) into SU(3). (In other words, in the present model, 
the flavor symmetry U(3) is explicitly broken from the begining by the Z2 symmetry. ) Under 
the requirement of the Z2 invariance, i.e. the invariance under the transformation 

(</.(«),</,«) ^ (-</.(«),+</.«), (2.5) 

the terms Tr(0(^)0(^V*'^'')) i-c- (~0^ + (l/3)</>^)^jj -\ , are forbidden. Thus, the superpotential 

(2.3) with the Z2 invariance leads to 
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The form (2.6) is just identical with (1.8) except for the "• • •" terms. As we show below, the 
potential (2.6) can give the desirable VEV relation (1.5) together with (</)]) = (i 7^ j). 

Prom the superpotential (2.6) with the Z2 invariance, we obtain the VEV relation (1.5) as 
follows: Prom the condition 

dW 



5(0(8) )■ 



= m 
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(2.7) 



we obtain 

m + X(j)a = 0, 

for {(t>^^^)l 7^ 0. By eliminating m from Eq.(2.8) and the condition 



(2.8) 



= 0, (2.9) 
we obtain the relation 

^l = Mcp('U^'^) = <t>l + ^l + ---, (2.10) 

where "• • • " denotes the contributions of 3 and 3' of S4. 

The result (2.10) is still not our goal, because the relation contains the VEVs of the 3 and 
3' of S4. So far, we have not discussed the splitting among the S4 multiplets. Now, we bring a 
soft symmetry breaking of SU(3) into S4 with an infinitesimal parameter e into the mass term 
of W{(f>) as 

TV(<^(«)0(«)) ^ <^,<^, + 4>r,4>r, + (!+£) J2{^^'^)i{4>^'% (2.11) 

by hand. (At present, we do not refer the origin of the symmetry breaking. The SU(3) flavor 
symmetry is explicitly (not spontaneously) broken with the order of e.) Recall that when we 
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obtain the relation (2.8), we have assumed {(f)^^^)l 7^ 0. Now, the conditions (2.7) are modified 
into the folloing conditons: 




Thefore, we must take either (^^^^)*- = (i 7^ j) or = (a = tt, rj) for £ 7^ 0. When we choose 
the solution 




we can obtain the desiarable relation (1.5). (However, it is possible that we can also take another 



is merely one of the possible solutions.) 

Thus, we have obtained not only the desirable VEV relation (1.5), but also the results 
(2.14). It should be worthwhile noticing that if we have assume the superpotential (2.3) without 
requiring the Z2 invariance, we could obtain neither (1.5) nor (2.14). 

3 Superpotential with symmetry breaking 

Since we know that the three masses in any sectors of quarks and leptons are completely 
different among them, we must consider that any flavor symmetry which we introduced should 
finally be broken completely. Although the superpotential (1.8) can give the VEV relation (1.5), 
it cannot fix the ratio Vtt/wjj. In order to fix the ratio v-j^/vrj, we consider the existence of an 
S4 symmetry breaking term WgB- Then, the problem is whether the VEV relation (1.5) which 
has been obtained under the S4 symmetry which is embedded into SU(3) is spoiled or not by 
introducing such a symmetry breaking, because, usually, a relation which we have derived under 
an exact symmetry is only approximately satisfied under the symmetry breaking. In the present 
section, we will demonstrate that such a symmetry breaking term without spoiling the relation 
(1.5) is indeed possible. 

We consider that the S4 invariant superpotential (1.8) is softly broken. Since we want 
v-^jvf^ 7^ 1, the breaking should appear in the doublet part of S4. In order to express the S4 
symmetry breaking term explicitly, we define the following symmetry breaking parameters B^^^ 
and with 3x3 matrix forms. 



solution with (p-,^ = = and 



(</.(8)) 



I*- 7^ 0. The VEV solutions are not unique. The result (1.5) 




(3.1) 



which behave as if those were octet and singlet of SU(3), respectively, where 



6^ = V2sin/3, 67r = V2cos/3, 6^ = 1, (3.2) 

and the factor \/2 in Eq.(3.2) has been chosen as 6^ + 6^ = 2 compared with 6^ = 1. Then, we 
can express the symmetry breaking term as the form 



WsB = ^em ['IV(S(«),^(«)0(«)) + TV(5(i)0(i)</.(i)) 



= [-2<f>^(f)r, COS P- {4>l-(f>l) sin P + (f>l] , (3.3) 

where the factor V3/2 has been chosen as the coefficients in the expression (3.3) correspond to 
those in the unbroken form (1.8). Although the term T}r{B^^') <f>^^') (f)^^^) = </>^/\/3 in (3.3) does 
not break the S4 symmetry, it has been added by hand in order that the term Wsb (in other 
words, the parameter e) does not affect the VEV relation (1.5). 

As the result, we can write the superpotential including the symmetry breaking term as 
follows: 

W = ^m {(t>l + + (1 + e)(^l - e [2cj)^(Pn cos/3 + {cj)l - (P^) sin/3] } + ^A^,, (^4>l + + ^0^) • 



Since 
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dW 

= [m + A</>o- — em sin /3] c/;^ — em(/)ri cos /3, (3-5) 
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dW 

— — = [m + + em sin /?] 0^ — emcp^^ cos /?, (3-6) 
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m(l + £)</., + -A((/>2 +02 + ^2)^ (3^7) 



the minimizing conditions of the potential leads to the relations 

v'^ — 

iim(3=— ^, (3.8) 

2VT,Vr^ 

vl + vl = vl, (3.9) 

m(l +e) + At;„ = 0. (3.10) 

Note that the derivation of the relation (3.8) is independent of the explicit values of m, A and 
e, and the derivation of the relation (3.9) is independent of the explicit values of m, A, e and 
/3. Thus, we can fix the value of Vt^Iv^ by the parameter /? in Wsb without spoiling the VEV 
relation (1.5) [(3.9)]. Also note that the limit m^. corresponds to the limit v^^ — > —v^/^/2 
(i.e. v"^ = u^), so that the limit m^ corresponds to /? ^ 0. 

When we define the parameters Zi = \Jm(.il ^m^ + m^^ + m,-, from the observed values [1] 
of the charged lepton masses, we obtain the numerical values z\ = 0.016473, Z2 = 0.236869 and 
Z3 = 0.971402, so that, for the VEVs of 4>a defined by Eq.(1.3) [(2.2)], we obtain z^ = 0.519393, 
-0.479824 and z^ = I/V2 = 0.707106. Therefore, we can estimate the value of (3 as 
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follows: 
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2 ^2 



smp= '^ ^ = 4z^ - 1 = -0.079078, /? = -4.5355°, (3.11) 



where we have chosen the phase convention of (3 as cos /3 = —2zTrZr^/{z'^ + ^V) > 0- 



Table 1 SU(3) and S4 assignments of the fields 
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From the point of view of the prameter physics, the introducing the symmetry breaking 
term (3.3) is merely replacing the parameter Utt/wjj by another parameter (3. What is important 
is that we can indeed introduce a symmetry breaking term without spoiling the relation (1.5). 

4 Effective Hamiltonian 

If we regard the scalars 4>u and 4>d as 'S\J(2)l doublets, such a model with multi-Higgs 
doublets causes a flavor changing neutral current (FCNC) problem. Therefore, we must consider 
that the fields (j)u and 0^ are SU(2)i singlets. In the present paper, we assume a Froggatt-Nielsen 
[12] type model 

H = y<^^LHL — — -eR + yjLHL — -VR + vrvr^rVr, (4.1) 

where In are SU(2)l doublet leptons liL = {i^iLj^iL), and H'^ are conventional SU(2)l 
doublet Higgs scalars, cpj {f = u,d), ^ and x are S\J{2)l singlet scalars, and A is a scale of the 
effective theory. We consider that {(pf)/A, (^)/A and {x)/^ are of the order of 1. The scalar ^r 
has been introduced in order to generate the Majorana mass Mr of the right-handed neutrino 
vr. As we note later, in the present model, the right-handed neutrinos ur = (i^j^^ + i^jj ^)/^/2 
are singlets of the SU(3) flavor. The role of = (^^^-^ + ^^~-*)/\/2 and x will be explained 
later. In order to understand the appearance of the combinations Hf^cpd^Pdi and Hf^cftuX, we 
assume two Z3 symmetries (Z3 and Z3 in Table 1). Those quantum number assignments are 
given in Table 1. However, even with those quantum numbers, we cannot distinguish the state 
from <pf<pf- For example, the interaction iLHd4>\^dR is possible in addition to the interaction 
^LHfi(t>d4'd^dR- Although we have started from an SUSY senario in the previous section, now, 
we have adopted an effective Hamiltonian which is not renormalizable. Therefore, in principle, 
the interaction (■LHd4>\^dR cannot be ruled out. For the moment, in order to forbid such an 
undesirable term, we assume that the fields which can appear in the effective Hamiltonian are 
confined to holomorphic ones. 



(a) Charged lepton sector 

Recall that we have already assumed the invariance of the superpotential under the Z2 
transformation (2.5) in order to drop the cubic part of the octet (f>^^\ Therefore, the term 



means (f)^^^(f)^^^ + under the Z2 invariancc. However, in order to give md oc 



what we want is not (jj^^^cp^^'^ + but ^(s)^^^) + ^(i)(/>(i) + ^(^V^^^ + ^^^V^^^- In order to 

evade this problem, we introduce additional fields ^^^^ and whose Z2 parity are +1 and —1, 
respectively. The effective interactions in the charged lepton sector are given by 

H!^^ = ^e-U0,)K0,),^(eW + &^)eRk, (4.2) 

where we have dropped the Higgs scalar Hf since we discuss flavor structure only. The expression 
(4.2) becomes 

= ^e-.[(<^f </>f + Vi^^)C(+^ + (tf Vi^^ + <t>'lUT)&W (4.3) 

Since we have assumed that ^^"^^ and appear symmetrically in the theory, we also assume 

{^i+)) = {^i-))^v^. (4.4) 
Then, we obtain the effective Hamiltonian for the charged leptons 



He 



= '-^T.^U{<t>'r')^'en. (4.5) 



where = {Hf). Since the fields {(t>d)\ are defined by Eq.(2.2), we can obtain the charged 
lepton mass relation (1.1) from the VEV relation (1.6). 

However, the present mechanism to obtain rngj oc is somewhat artificial. The present 
mechanism will be improved in the future model. [Of course, there is a possibility that the 
superpotential (2.3) must exactly be invariance under the Z2 symmetry, but the effective Hamil- 
tonian (4.1) does not need to be invariance under the Z2 symmetry. Then, we can consider a 
model without .] 

(b) Neutrino sector 

In the present model, the right-handed neutrinos v^'^^ are singlets of SU(3). Therefore, in 
the neutrino seesaw mass matrix Mi, = m'£M^^{m'[^)^ , Mr is a 1 x 1 matrix and is a 3 x 1 
matrix. In order to compensate for the absence of the conventional triplet neutrinos vr, a new 
scalar % which is a triplet of SU(3) has been introduced. The neutrino Dirac mass terms are 
given by the following effective Hamiltonian 

<L = y.^^i(('/'«)^)(x,)(4^^ + 4"^ (4-6) 

where Vu = {HJ^). It is likely that the scalar potential V{x) the SU(3) triplet x has a specific 
VEV solution 

(Xi) = (X2) = (X3) = vx- (4.7) 



When we assume the VEVs (4.7), we obtain 
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where Va = {(^ua) (« = ^r, 77, a) (for convenience, we have dropped the index u). Therefore, we 
obtain the effective neutrino mass matrix on the (r/, ex, tt) basis. 
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where = ($_r), and we have dropped the common factors [yyVuV^j ^f^b^)"^ . By the 

way, the ratio Vt^/v.^ cannot be determined from the potential (2.6), and the ratio is determined 
by a soft S4 symmetry breaking term Wsb which has been discussed in the previous section. 
We can choose a solution = Q \n the superpotential W{4)u) by adjusting the parameter j3 
in WsBi differently from the case of W{(j)d)- Then, the neutrino mass matrix (4.9) becomes a 
diagonal form = (1/M^" Viag(?;2, 0, 0) + (l/M^"^Viag(0, f^, 0). Since the mass matrix My 
on the (i^i, z/2, u-^) = (z^e> ^'/x, i^t) basis is given by 

My = UtbMI'"'^^ U^b = UTBDyU^B, (4.10) 

we can obtain the tribimaximal mixing 

Uy = Utb, (4.11) 

and the neutrino masses 

ruyi = kv"^, mu2 = kvl, 171^3 = 0, (4.12) 

for the case of M^j^^ = M^j ^ = AIr, where k = (yyVuV^)'^ /2MfiA'^ and (fr;, fa, i^tt) has been 
renamed (i^i, 1^2, t'a) according to the conventional naming. 

However, since we have taken v^^ = 0, the value of v^^ satisfies v'^ = from the relation 
(1.5), so that the result (4.12) gives rriyi = 171^2- The observed value [H] Am^^;^^ is small, but 
it is not zero. Therefore, we must consider a small deviation between the first and second terms 
in (4.9) (i.e. M^j^^ ^ M^ ^). Since the value /M^^ is free in the present model, we cannot 
predict an explicit value of the ratio Am^^j^^/Am^^^. 

Since the present model gives an inverse hierarchy of the neutrino masses, the predicted 
effective electron neutrino mass 



E 



U^iTJly 



\m 
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^•^■^-0.40 ^ 



10"^ eV, 



(4.13) 



where we have used the value [15] Am^tm = 2.74]1^q;26 x 10 ^ eV^. This value (4.13) is sufficiently 
sensitive to the next generation experiments of the neutrinoless double beta decay. 

5 Summary 

In conclusion, on the basis of the S4 symmetry which is embedded into SU(3), we have 
investigated a lepton mass model with the effectuve Hamiltonian of the Froggatt-Nielsen type 
(4.1). We have assumed that the singlet and doublet of S4 originate in the singlet and octet 
of SU(3), and we have obtained the VEV relation (1.5). In the derivation of the VEV relation 
(1.5), the essential assumptions for the superpotential W{(j)f) are the following two: (i) the scalar 
fields (j)f always appear in terms of the nonet form (2.1) of U(3); (ii) the superpotential W{(j)f) is 
invariant under the Z2 transformation (2.5). Then, we have obtained not only the VEV relation 
(1.5), but also ((</'(^))-) = (i / j) for the other components of ^^^^ (i.e. (3) = (3') = 0). 

In the charged lepton secter, in order to give rriei oc {{4'd)i)'^j we have assumed new scalars 
Although it has been reuired to compensate for the Z2 invariance, the model seems to 
leave the door open to further improvement. 

For the neutrino sector, we have obtained the tribimaximal mixing (1.2) by introducing an 
SU(3) triplet scalar x and the two SU(3) singlet right-handed neutrinos u^^^ in addition to the 
nonet scalar In the present model, the right-handed neutrinos v^^^ are singlets of SU(3), 
the Majorana neutrino mass matrices M^'^ have no flavor structure. For the neutrino mass 
spectrum, since the model gives m,yi = my2 in the limit of m{j+^ = \ 

we must consider a 

small deviation Afjj''^ / Since the value of M^j^V^/T^ 

is a free parameter in the present 
model, we cannot predict the value Am^^j^^/Am^^^ at present, although the smallness of the 
ratio Am^^i^,^ / Am'^fjj^ can be understood. Since the present model gives an inverse hierarchy of 
the neutrino masses, we can predict the effective electron neutrino mass {m^^) ~ 0.05 eV, which 
is sufficiently sensitive to the next generation experiments of the neutrinoless double beta decay. 

The present model seems to provide suggestive hints on seeking for a model which leads to 
the tribimaximal mixing (1.2) and the charged lepton mass relation (1.1), although the model 
has still many points which should be improved. At the same time, the model will provide a 
clue to the quark mass matrix model from a point of unified view of the quarks and leptons. 
The extension of the present model to the quark mass matrix model will be given elsewhere. 

Acknowledgements 

The author would like to thank E. Takasugi, T. Fukuyama and H. Fusaoka for helpful 
conversations. Especially, the author is much indebted to N. Haba for his valuable contribution 
to the improvement on the previous version. He also indebted to H. Fusaoka for the phase 
convention of S4. This work is supported by the Grant-in-Aid for Scientific Research, Ministry 
of Education, Science and Culture, Japan (No. 18540284). 



References 



[1] Particle Data Group, J. Phys. G: Nucl. Phys. 33 (2006) 1. 



[2] Y. Koide, Lett. Nuovo Cimento 34 (1982) 201; Phys. Lett. B120 (1983) 161; Phys. Rev. 
D28 (1983) 252. 



[3] Y. Koide, Mod. Phys. Lett. A5 (1990) 2319. 

[4] S. Pakvasa and H. Sugawara, Phys. Lett. B82 (1979) 105; Y. Yamanaka, H. Sugawara and 
S. Pakvasa, Phys. Rev. D25 (1982) 1895; P. F. Harrison, D. H. Perkins and W. G. Scott, 
Phys. Lett. B458 (1999) 79; Phys. Lett. B530 (2002) 167; Z. Z. Xing, Phys. Lett. B533 

(2002) 85; P. F. Harrison and W. G. Scott, Phys. Lett. B535 (2003) 163; Phys. Lett. B557 

(2003) 76; E. Ma, Phys. Rev. Lett. 90 (2003) 221802; C. I. Low and R. R. Volkas, Phys. Rev. 
D68 (2003) 033007; X.-G. He and A. Zee, Phys. Lett. B560 (2003) 87; R. N. Mohapatra, 
S. Nasri and H. B. Yu, Phys. Lett. B639 (2006) 318. 



[8. 
[9 

[lo; 
[11 
[12; 

[13 
[14 



Y. Koide and H. Fusaoka, Z. Phys. C71 (1996) 459. 
Y. Koide and M. Tanimoto, Z. Phys. C72 (1996) 333. 

S. Pakvasa and H. Sugawara, Phys. Lett. B73 (1978) 61; H. Harari, H. Haut and J. Weyers, 
Phys. Lett. B78 (1978) 459; J.-E. Erere, Phys. Lett. B80 (1978) 369; E. Derman, Phys. Rev. 
D19 (1979) 317; D. Wyler, Phys. Rev. D19 (1979) 330. 

Y. Koide, Phys. Rev. D60 (1999) 077301. 

Y. Koide, Phys. Rev. D73 (2006) 057901. Also, see Y. Koide, Phys. Rev. D60 (1999) 
077301. 

Y. Koide, arXiv: |hep-ph/0701018l 
E. Ma, arXive: |hep-ph/ 0612022[ 

C. Froggatt and H. B. Nielsen, Nucl. Phys. B147 (1979) 277. 

C. Hagedorn, M. Linder and R. N. Mohapatra, JHEP, 0606 (2006) , 042 |hep-ph/0602244l . 

B. Aharmim et al. SNO Collaboration, Phys. Rev. C72 (2005) 055502. T. Araki et al. 
KamLAND Collaboration, Phys. Rev. Lett. 94 (2005) 081801. 



[15] D. G. Michael et al. MINOS Collaboration, Phys. Rev. Lett. 97 (2006) 191801. Also see, 
J. Hosaka et al. the Super-Kamiokande Collaboration, arXiv: hep-ex/06040TT| 



